arXiv:l502.00643v2 [astro-ph.CO] 1 Sep 2015 


TO BE SUBMITTED TO ASTROPHYS. J.— DRAFT VERSION SEPTEMBER 2, 2015 
Preprint typeset using style emulateapj v. 5/2/11 


Bicep2 / Keck Array V: MEASUREMENTS OE B-mode POLARIZATION AT DEGREE ANGULAR SCALES AND 

150 GHZ BY THE Keck Array 

Keck Array AND BICEP2 COLLABORATIONS: P. A. R. Ade', Z. Ahmed^, R. W. Aikin^, K. D. Alexander"^, D. Barkats^, 

S. J. Benton'’, C. A. Bischoef"', J. J. Bock^’^, J. A. BREVIK^ I. Buder"^, E. Bullock*, V. Buza"*, J. Connors'*, B. P. Crill* ’, 
C. D. Dowell’, C. Dvorkin”*, L. Duband^, J. P. FilippinE ''*, S. Fliescher", S. R. Golwala*, M. Halpern*’, 

M. Hasseleield*’, S. R. Hildebrandt*’’, G. C. Hilton'*, V. V. Hristov*, H. Hui*, K. D. Irwin* '*’''*, K. S. Karkare'*, 

J. P. Kaufman'*, B. G. Keating'*, S. Kefeli*, S. A. Kernasovskiy*’’^, J. M. Kovac**, C. L. Kuo*’*'*, E. M. Leitch'®’*’, 

M. Lueker*, P. Mason*, K. G. Megerian’, C. B. Netterfield* '*, H. T. Nguyen’, R. O’Brient’, R. W. Ogburn IV* *'*, 

A. Orlando'*, C. Pryke"’*, C. D. Reintsema**, S. Richter**, R. Schwarz", C. D. Sheehy'* *’, Z. K. Staniszewski* ’, 

R. V. Sudiwala', G. P. Teply*, K. L. Thompson*, J. E. Tolan*, A. D. Turner’, A. G. Vieregg'* '’, A. C. Weber’, 

J. Willmert", C. L. Wong'*, and K. W. Yoon*-''* 
to be submitted to Astrophys. J. — draft version September 2, 2015 

ABSTRACT 

The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the 
Bicep2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which 
the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing 
field as Bicep2. We again find an excess of B-mode power over the lensed-ACDM expectation of > 5(t in the 
range 30 < f < 150 and confirm that this is not due to systematics using jackknife tests and simulations based 
on detailed calibration measurements. In map difference and spectral difference tests these new data are shown 
to be consistent with Bicep2. Finally, we combine the maps from the two experiments to produce final Q and 
U maps which have a depth of 57 nKdeg (3.4 pKarcmin) over an effective area of 400 deg^ for an equivalent 
survey weight of 250,000 /iK“^. The final BB band powers have noise uncertainty a factor of 2.3 times better 
than the previous results, and a significance of detection of excess power of > 6cr. 

Subject headings: cosmic background radiation — cosmology: observations — gravitational waves — infla¬ 
tion — polarization 


1. INTRODUCTION 

Precision polarimetry of the cosmic microwave background 
(CMB) has become a mainstay of observational cosmology. 
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The ACDM model predicts a polarization of the CMB at 
the level of a few /xK, with a characteristic B-mode pat¬ 
tern. The EE power spectrum has been detected over a 
wide range of angular scales by many experiments, includ¬ 
ing Dasi (Kovac et al. 2002), CAPMAP (Barkats et al. 2005; 
Bischoff et al. 2008), CBI (Readhead et al. 2004; Sievers et al. 
2007), Boomerang03 (Montroy et al. 2006), WMAP (Page 
et al. 2007; Bennett et al. 2013), MAXIPOL (Wu et al. 
2007), quad (Pryke et al. 2009; Brown et al. 2009), Bl- 
CEPl (Chiang et al. 2010; BlCEPl Collaboration et al. 2014), 
Quiet (Quiet Collaboration et al. 2011, 2012), Polar- 
bear (POLARBEAR Collaboration et al. 2014a), Bicep2 (Bi- 
CEP2 Collaboration I 2014), ACTPol (Naess et al. 2014), 
SPTpol (Crites et al. 2014), and Planck (Planck Collabora¬ 
tion results I 2014). These measurements have been in broad 
agreement with theoretical expectations and other cosmolog¬ 
ical data sets. Improved EE power spectrum data are impor¬ 
tant because they may eventually constrain the ACDM model 
parameters better than cosmic variance limited CMB temper¬ 
ature data (Rocha et al. 2004; Galli et al. 2014). 

Of greater interest is the B-mode component of the polariza¬ 
tion pattern. Though the EE power spectrum is higher, the BB 
power spectrum is more sensitive to new physics because the 
linear density perturbations at the surface of last scattering, 
which are the main source ofTT and EE power, cannot gener¬ 
ate B-mode power. On small angular scales, BB power instead 
arises from the gravitational lensing of B-mode power by the 
large scale structure of the universe (Zaldarriaga & Seljak 
1998). The lensing BB power thus cleanly traces the growth of 
structure, complementary to other methods, providing infor¬ 
mation about possible extensions to ACDM such as neutrino 
mass or a nontrivial dark energy equation of state. Measure- 
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ments by SPTpol (Hanson et al. 2013) POLARBEAR (POLAR- 
BEAR Collaboration et al. 2014b,c,a), and Bicep2 (Bicep 2 
Collaboration I 2014) have provided the first evidence of BB 
power from gravitational lensing. 

On large angular scales, lensing contributes only a small 
amount of BB power. However, inflationary gravita¬ 
tional waves (IGW) may be a source of BB power on 
these scales (Polnarev 1985; Seljak & Zaldarriaga 1997; 
Kamionkowski et al. 1997; Seljak 1997). The recent detec¬ 
tion by Bicep2 of B-mode power on degree angular scales 
in excess above the lensing expectation is especially excit¬ 
ing because it could be evidence of primordial gravitational 
waves and cosmic inflation (Bicep2 Collaboration I 2014). 
The contribution of foregrounds to the observed Bicep2 sig¬ 
nal is uncertain, and preliminary data from Planck have sug¬ 
gested that polarized dust in the Bicep2 field may be brighter 
than models had predicted (Planck Collaboration Int. XXX 
2014). Regardless, a confirmation of the Bicep2 signal, 
whether cosmological or galactic in origin, is a top priority 
of observational cosmology today (Dodelson 2014; Caligiuri 
& Kosowsky 2014). 

The Keck Array telescope is a microwave polarimeter at the 
South Pole designed to follow up the Bicep 2 observations. 
The Keck Array quickly deployed a large number of detec¬ 
tors at 150 GHz by installing five receivers of a design very 
similar to Bicep 2 with minimal changes. All five receivers 
were installed in time for the 2012 observing season and con¬ 
tinued, with upgrades, to observe at 150 GHz through the end 
of 2013. The modular, multi-receiver structure of the Keck 
Array also allows individual receivers to be tuned to different 
frequencies. Two of the Keck Array receivers began observing 
at 95 GHz in 2014, which will help to discriminate the signal 
from foregrounds. 

In this paper, we present the results of the 150 GHz ob¬ 
servations by the Keck Array of the Bicep2 field during the 
2012 and 2013 seasons. We begin with sections describing the 
Keck Array instrument, calibrations, and analysis methods. 
We proceed with the maps and angular power spectra obtained 
from this data set and perform extensive internal consistency 
checks. The Keck Array confirms the Bicep2 B-mode signal 
at > 5 ct. We then proceed to test for consistency between the 
Bicep2 and Keck Array data, and finally combine the two sets 
of maps to a final result. 

This paper is the latest in a series of publications by the Bl- 
CEP2 and Keck Array collaborations. Bicep2 Collaboration I 
(2014) (hereafter the Bicep2 Results Paper) is directly anal¬ 
ogous to this paper. Bicep2 Collaboration II (2014) (here¬ 
after the Bicep2 Instrument Paper) presented the full details 
of the Bicep2 instrument—differences are summarized in §2 
of this paper. Bicep2 Collaboration III (2015) (hereafter the 
Systematics Paper) presents a detailed analysis of instrumen¬ 
tal systematics, which are treated for the Keck Array 2012- 
13 data in §7 of the current paper. Two additional papers, 
Bicep2 / Keck Array Collaborations IV (2015) (the Beams 
Paper) and Bicep2, Keck Array, and SPIDER Collaborations 
(2015) (the Detectors Paper), describe the beam characteriza¬ 
tion and the detectors for both Bicep2 and the Keck Array. 

2. THE Keck Array INSTRUMENT 

The Keck Array instrument shares much of its design with 
Bicep2, details of which are presented in the Bicep2 Instru¬ 
ment Paper. In this section, we describe the main features 
common to both instruments and the substantive changes and 
upgrades unique to the Keck Array. Figure 1 shows the re¬ 


ceiver design for the Keck Array. 

2.1. Cryostat and cryogenic system 

The Keck Array comprises five independent 
cryostats (Sheehy et al. 2010) built by Atlas Technolo- 
gies^°. Inside each cryostat is a closed-cycle, three-stage 
("^He/^He/^He) sorption refrigerator (Duband (& Collaudin 
1999) that cools the focal plane unit (FPU) to approximately 
270 mK. Other optical elements are held at cryogenic 
temperatures to minimize the thermal load on the FPUs. 

The main difference between Keck Array and Bicep 2 is 
the bulk refrigeration system. While Bicep 2 used a bath of 
liquid helium, the Keck Array uses a set of Cryomech^' PT- 
410 pulse tube refrigerators. Each Keck Array cryostat has 
its own pulse tube refrigerator aligned along the optical axis. 
The helium gas is pulsed at a common frequency of 1.2 Hz, 
and the pressure in each system is optimized to achieve the 
lowest base temperature. After optimization, the pulse tubes’ 
copper mounting surfaces typically reach 40 K and 3 K with 
comparable performance in all five cryostats. These surfaces 
are thermally connected to the telescope insert by stacks of 
ultra high purity aluminum foil. 

2.2. Optics 

The Keck Array optics use an on-axis, refractive de¬ 
sign which was originally demonstrated in the BlCEPl tele¬ 
scope (Takahashi et al. 2010). The entire optics chain is es¬ 
sentially unchanged from Bicep2 (Aikin et al. 2010). The 
two-lens design was chosen to accommodate the flat telecen- 
tric focal plane, with the image of the primary at infinity as 
viewed from the focal plane. 

The lenses are made from high density polyethylene 
(HOPE) and cooled to 4 K. In order to sufficiently reduce the 
infrared loading on the cooling stages, there is a 3 mm nylon 
filter and two polytetrafluoroethylene (PTEE) filters of thick¬ 
ness 12.7 mm and 34.3 mm in the optics path, cooled to an 
intermediate stage of 50 K. A second nylon filter of 5.2 mm 
is placed between the objective and eyepiece lens, heat sunk 
to 4 K. A metal mesh low-pass filter (Ade et al. 2006) with a 
cutoff of 8.3 cm“' (225 GHz) is placed above the final nylon 
filter to prevent any stray radiation that was not absorbed by 
the plastic filters from thermalizing in the detectors. 

All surfaces surrounding the optical path are blackened with 
Eccosorb HR 10^^ cut in half and epoxied with Stycast 2850 
loaded with carbon. The Stycast covers the HRIO to prevent 
particulate shedding during cryogenic cycling. The lining is 
designed such that stray light terminates on cold surfaces. 
Later configurations installed after the 2013 observing season 
included baffling on the inside of the telescope tube to further 
reduce reflections. 

The aperture stop is similarly made from a ring of 1.9 cm 
thick Eccosorb AN-74, beveled at 40° with an inner diameter 
of 26.4 cm. Approximately 20% of the total throughput is 
absorbed by the aperture stop. The aperture stop is placed on 
the lower surface of the objective lens. 

2.3. Focal plane unit 

The detectors, developed at Jet Propulsion Laboratory 
(JPL) for joint use in the Bicep2, SPIDER and Keck Array 

^*^http: //www. atlasuhv. com/ 

^%ttp://www.cryomech.com/ 

^^Emerson & Cuming Microwave Products, Randolph, MA 02368. 
Phone: 781-961-9600. Web: http: //eccosorb.com/ 
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Fig. 1.— Individual receiver of the Keck Array. Each receiver is cryogenic, with a pulse tube refrigerator cooling the optics to 4 K and a three-stage sorption 
refrigerator cooling the focal plane to 270 mK. The Keck Array consists of five identical receivers on a single telescope mount at the South Pole. 


experiments, consist of planar antennas, each an array of slot 
sub-antennas combined in phase, and feeding into dual series 
Ti and A1 Transition-Edge Superconducting (TES) bolome¬ 
ters. A detailed description is given in the Detectors Paper. 
Each pixel consists of two interleaved phased antenna arrays 
with orthogonal polarization directions. The signal is band¬ 
pass filtered by a lithographed filter in stripline after the an¬ 
tenna and terminated on a thermally isolated island which also 
contains the series TES. A single Si tile contains 64 detector 
pairs and a focal plane unit has 4 tiles. In each focal plane, 
8 detector pairs are left “dark”. Dark detectors consist of the 
complete TES island structure, but are not connected to their 
corresponding antennas. 

The focal plane units in the Keck Array were slightly mod¬ 
ified versus those described in the Bicep2 Instrument Paper. 
The spacing between the tiles was increased in order to reduce 
the electromagnetic coupling between the pixels near the edge 
of the tile and the copper plate. The feed network of the anten¬ 
nas was also redesigned to reduce the dipole beam mismatch 
between the polarized pairs, significantly improving the per¬ 
formance versus Bicep2. 

The TES detectors are voltage biased and the current is in¬ 
ductively coupled to time-domain multiplexing SQUIDs (de 
Korte et al. 2003). The Keck Array uses the NIST developed 
MUX09s, which have a gradiometric design that reduces the 
sensitivity of the SQUIDs to uniform magnetic fields by three 
orders of magnitude in comparison to the MUX07a design 
used in Bicep2 (Stiehl et al. 2011). 

After the deployment of Bicep2, we discovered that the 
aliased noise from the multiplexing system was affecting the 
overall sensitivity of the instrument. One way to mitigate the 


aliased noise is to increase the Nyquist inductors that limit 
the bandwidth of the detectors. The choice of inductance is 
balanced with the need for the L/R time constant to be fast 
enough for the detectors to be in stable negative electrother¬ 
mal feedback. The first focal plane produced for the Keck 
Array has Nyquist chips with an inductance of 1.35 pH con¬ 
sistent with Bicep2, and all subsequent focal planes have an 
increased inductance of 2 pH. This limits the bandwidth to 
< 18 kHz. 

2.4. Readout 

The configuration of the room temperature electronics that 
interface with Keck Array is similar to that described in 
the Bicep 2 Instrument Paper. A Multi-Channel Electronics 
(MCE) crate provided by the University of British Columbia 
mounts directly to the outside of each cryostat to interface 
with the SQUIDs, and supply the detector bias (Battistelli 
et al. 2008). The Keck Array MCE crates use lower power 
SQUID series array readout cards compared to earlier designs 
in order to stay compatible with a development program to 
improve operability on balloon-borne telescopes. 

Housekeeping thermometry is read through “backpacks” 
attached to the cryostats similar to Bicep 2, and the signals 
are collated and digitized in a common BLASTbus2 crate pro¬ 
vided by the University of Toronto (Benton et al. 2014). Both 
the housekeeping and detector electronics are connected to a 
set of Linux-based computers and recorded to disk using the 
control software gcp at a sample rate of 20 Hz (Story et al. 
2012 ). 

2.5. Mount 
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TABLE 1 

The Keck Array DETECTOR PARAMETERS 


Detector Parameter 

Median Value 

Optical efficiency, r] 

24% 

Band center, (i/) 

151 GHz 

Spectral bandwidth, Au 

42 GHz 

Normal resistance, 

62 ml2 

Operating resistance, R^p 

0.68 P/v 

Saturation power, Psat 

9.9 pW 

Optical loading, Popt 

3.1 pW 

Thermal conductance, Gc 

90 pW/K 

Transition temperature, Tc 

520 mK 

Thermal conductance exponent, (3 

2.5 


All five Keck Array receivers are attached to a common 
telescope mount located at the Martin A. Pomerantz Obser¬ 
vatory (MAPO) at the Amundsen-Scott South Pole Station. 
This mount was previously used for the Dasi (Leitch et al. 
2002) and QUaD (Hinderks et al. 2009) experiments. A new 
front end “drum” for the Keck Array cryostats was installed 
in 2010. The platform was leveled at that time to account for 
gradual shifts of the building on the snow relative to the hori¬ 
zon. 

The mount has three axes: elevation, azimuth, and bore- 
sight. The rotation around the boresight is referred to as “deck 
rotation” and allows for cancellation of some systematic ef¬ 
fects and/or tests for their presence. 

3. CHARACTERIZATION 

The Keck Array has been extensively characterized in labo¬ 
ratory tests and with in situ calibration measurements. The 
characterization program was very similar to the one de¬ 
scribed in Sections 10 and 11 of the Bicep 2 Instrument Pa¬ 
per. This section summarizes these measurements, focusing 
particularly on detector properties that have been reoptimized 
since the fabrication of the Bicep 2 detectors. The spectral 
band, optical efficiency, and bolometer thermal conductance 
have a strong effect on the ultimate sensitivity of the instru¬ 
ment, and are tuned to minimize noise while allowing stable 
operation under typical South Pole atmospheric loading con¬ 
ditions. These detector properties, summarized in Table 1, are 
described in Sections 3.1, 3.2, and 3.3. We have also exten¬ 
sively measured the far-held beams with in situ observations 
of a mast-mounted source. The beam mapping measurement 
and its results are summarized here in Section 3.4 and de¬ 
scribed more fully in the Beams Paper. 


3.1. Spectral response 

The frequency response of the antennas and lumped el¬ 
ement filters was tuned to give a fractional bandwidth of 
25%. The ~ 150 GHz observing band is bracketed by the 
118.8 GHz oxygen line on the low side and the 183.3 GHz 
water line on the high side. 

The frequency response S{v) of the 150 GHz detectors was 
characterized using a Martin-Puplett Fourier transform spec¬ 
trometer (FTS) (Karkare et al. 2014). From these spectra, the 
band center and the bandwidth are calculated. The band cen¬ 
ter is defined to be 


/ vI{v)S{i')dv 
f /(iy)S(iy)diy 


(1) 


and the bandwidth 

_ (f I(u)S(i^)di^f 
f K(iy)S^(i^)di' ’ 


(2) 


where /(i/) is the source spectrum relative to a Rayleigh-Jeans 
spectrum. For a Rayleigh-Jeans spectrum, the 150 GHz detec¬ 
tors were measured to have a band center of 151.5±1.9 GHz 
and a bandwidth of 41.8±1.4 GHz. The effective band center 
shifts to 150.6, 152.8 and 148.6 GHz respectively for a source 
spectrum of CMB, dust, and synchrotron radiation using the 
current best-fit models (Planck Collaboration Int. XXX 2014; 
Bennett et al. 2013) The standard deviations are dominated by 
variation from tile to tile, with smaller variation from detector 
to detector within a tile. 


3.2. Optical efficiency 

The end-to-end optical efficiency is defined as the fraction 
of incident light absorbed by the detectors. This is dependent 
on the losses in the optics, the antennas, and the bandpass 
filters. A higher optical efficiency increases the sensitivity of 
the detectors. Because it also increases the optical loading and 
photon noise, it must be taken into account when optimizing 
the thermal conductivity of the detector. 

The optical efficiency is measured in the lab using a beam¬ 
filling, microwave absorbing cone of AN-72. The power 
change on the detector for a source in the Rayleigh-Jeans limit 
(hi/ kT) is: 

Popt = kT pAi/ (3) 

where rj is the optical efficiency and Au is the bandwidth as 
defined in §3.1. The detector loading was measured with the 
cone at both room temperature and liquid nitrogen tempera¬ 
ture and converted into optical efficiency with the measured 
bandpass of 42 GHz. 

The median measured end-to-end optical efficiency of the 
Keck Array receivers was 24%. The optical efficiency of the 
early Keck Array detectors was lower than the 38% optical ef¬ 
ficiency of the Bicep2 detectors. Detector testing after initial 
Keck Array deployment suggested that the optical efficiency 
was being reduced by microscopic stress-induced cracks in 
the niobium microstrips connecting the antenna networks to 
the TES bolometers. In later generations of Keck Array de¬ 
tectors the Nb film stress was decreased from ^ 1000 MPa to 
< 300 MPa. The optical efficiency was observed to increase 
as described in more detail in the Detectors Paper. 

3.3. Thermal conductance 

As is discussed in the Bicep2 Instrument Paper, the detec¬ 
tor parameters can be tuned during fabrication in order to op¬ 
timize the noise performance. In particular, the thermal con¬ 
ductance can be tuned to minimize the phonon noise while 
maintaining a margin of safety ensuring operability under nor¬ 
mal loading conditions. 

The phonon noise is the thermal fluctuations from the sub¬ 
strate to the detector island through the SiN isolation legs. 
The noise-equivalent-power (NEP) is dependent on the ther¬ 
mal conductance G across the legs (see e.g. Irwin & Hilton 
(2005)) as 

NEP= ^AkeGT^F. (4) 

where E is a numerical factor describing the non-linearity of 
the thermal conductance between the substrate temperature 
and the bath temperature (typically 0.5 for these detectors). 
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The saturation power of the detectors is dependent on the 
thermal conductance as: 


■■GqTq- 


(r, /7b/+‘-i 

/ 3+1 


(5) 


where the exponent j3 is roughly 2.5 for these detectors. For 
the Keck Array the loading from the optics and the sky was 
modeled to be ~22 K„. The optical efficiency is used to con¬ 
vert the loading temperature to a power deposited on the de¬ 
tector. For the median optical efficiency described in §3.2, 
this corresponds to Popt =3.1 pW of loading under normal ob¬ 
serving conditions. Assuming a safety factor of 2, the optimal 
Gc is then 67 pW/K. 

The thermal conductance Gc was measured using detector 
load curves with the substrate held at different temperatures. 
This method used “dark” detectors that were purposefully dis¬ 
connected from their antennas to avoid the optical loading ef¬ 
fects. The detectors used in Bicep 2 had higher Gc, with two 
tiles centered at 100 pW/K and two centered at 140 pW/K. 
The tiles fabricated in later runs for the Keck Array had lower 
thermal conductances, with a median Gc of 90 pW/K. Several 
tiles had a much lower Gc of 30-50 pW/K, expected to give 
lower phonon noise but a smaller margin of safety against sat¬ 
uration. 

Finally, the margin of safety can be verified by measur¬ 
ing the electrical power Pj required to keep the detector in 
transition during standard observation. The standard observ¬ 
ing schedule includes load curves (bolometer current-voltage 
measurements) taken once per hour. These have been used 
retrospectively to assess the safety margin under actual atmo¬ 
spheric conditions. With the telescope pointed at 55° in ele¬ 
vation, the detectors were found to have a median margin of 
safety of 6.8 pW, corresponding to a safety factor of 3.2. 


TABLE 2 

The Keck Array OBSERVATION PHASES 


Phase 

LST Time 

Field 

Elevation [deg] 

Azimuth [deg] 

A 

Day 0 23:00 

Cryo service 



B 

Day 1 05:30 

CMB 

55.00-57.25 

120-300 

C 

Day 1 14:30 

CMB 

57.50-59.75 

-10-170 

D 

Day 1 23:00 

Galaxy 

55.00-56.50 

130-270 

e 

Day 2 05:30 

CMB 

57.50-59.75 

120-300 

F 

Day 2 14:30 

CMB 

55.00-57.25 

-10-170 


relative gain of the detectors. The scan rate in azimuth was 
2.8°/s. 

A group of ten scansets over successive azimuth ranges 
(and stepping in elevation) is called a “phase”. Table 2 shows 
the phases for the Keck Array. The elevation ranges were 
switched between phases after each full cycle of schedules. 
Since the briefest sub-kelvin hold time among the set of five 
helium sorption refrigerators was ~48 hours, the standard 
observing schedule consisted of four CMB phases and one 
galaxy phase between fridge cycles. 

As for Bicep2, the Keck Array mount allows for rotation 
of the whole apparatus around the line of sight—referred to 
as “deck rotation”. This rotation was performed between each 
two day schedule. For Bicep2 and the Keck Array in 2012, 
four deck angles were used: 68, 113, 248, and 293. These four 
angles provide coverage in Q and U and allow for cancellation 
of systematic effects whose sign reverses under 180° rotation. 
In 2013, the Keck Array started observing at eight deck an¬ 
gles: 23, 68, 113, 158, 203, 248, 293, and 338°. This allows 
for a more complete cancellation of beam systematics—see 
the Systematics Paper. 


3.4. Beams 

The beam shapes were measured in situ at the South Pole 
by scanning on a large thermal noise source mounted ~ 200 m 
away, in the optical far field. All receivers in the Keck Array 
have a beam width of 0.22 degrees, with very low levels of el- 
lipticity. As in Bicep2 the dominant differential beam imper¬ 
fection for the Keck Array is differential pointing. The beam 
mapping campaign, extracted beam parameters, and residual 
beam features are described in detail in the Beams Paper. In 
this paper we use the high-fidelity per-detector beam maps as 
a convolution kernel for simulations to place a limit on the 
false B-mode signal from beam imperfections. The simula¬ 
tions and results are described in §7. 

4. OBSERVATIONS AND DATA SET 

4.1. Observations 

The observation strategy of the Keck Array in the 2012 and 
2013 seasons was very similar to that used by Bicep2, as de¬ 
scribed in the Bicep2 Results and Instrument Papers. The 
same field as Bicep2 (and BlCEPl) was observed—a region 
centered at RA Oh, Dec. -57.5°. As viewed from the South 
Pole, the observing field remains at constant elevation and ro¬ 
tates in azimuth once per day. For fifty minutes periods the 
telescope scanned in azimuth at a fixed elevation, forming a 
“scanset” with 102 half-scans. Between scansets, the azimuth 
was updated by approximately 12.5° to account for the sky 
rotation, and stepped in elevation by 0.25°. Before and after 
each scanset, an elevation nod was performed to calibrate the 


4.2. Data selection 

As described in the Bicep2 Results and Instrument Papers, 
data cuts are applied in three distinct stages. A few cuts re¬ 
move half-scans from the scansets, while a larger number cut 
entire scansets from the final map. The final cut stage is the 
channel selection cut which is applied during the final coaddi¬ 
tion stage. These three stages provide the necessary flexibility 
and granularity. The cuts are summarized in Table 3 for the 
2012-2013 data set (Cf. Table 7 of the Bicep2 Instrument 
Paper). 

4.3. 2012-2013 data set and sensitivity 

The telescopes continuously took data through the South 
Pole winter. Each of 2012 and 2013 produced nearly 4500 
scansets of data for a total of 18 x 10^ s of data. After the data 
selection cuts, an average of 50% of the data remain, for a 
total of 15.6 X 10^ dets-s. This is shown graphically in the top 
panel of Figure 2 (Cf. Figure 23 of the Bicep2 Instrument 
Paper). 

The instantaneous sensitivity of the Keck Array instruments 
is measured using two methods: by taking the average of the 
time stream noise spectra between 0.1-1 Hz and by measur¬ 
ing the standard deviation of the noise-only maps weighted by 
the square root of integration time (Kernasovskiy et al. 2012). 
Both methods yielded a noise-equivalent-temperature (NET) 
of 11.5 for 2012 and 9.5 pK^.^s for 2013. The 

middle panel of Eigure 2 shows the instantaneous sensitivity 
calculated with the time stream based method for the 2012- 
2013 seasons. Using the same method as described in the 
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TABLE 3 

Data selection cuts EOR the Keck Array EOR 2012-2013 


Cut 

Total time [10® s] 

Integration [10^ det-s] 

Fraction cut [%] 

Before cuts 

18.3 

30.4 


Channel cuts 

18.3 

27.5 

9.5 

Synchronization 

18.1 

27.2 

1.2 

Deglitching 

17.8 

23.6 

11.8 

Passing channels 1 

17.6 

23.3 

0.58 

Elnod calibration 

17.1 

19.9 

11.2 

TES fractional resistance 

17.1 

19.6 

0.88 

Time stream skewness 

17.1 

17.8 

6.2 

Time stream variance 

17.0 

17.4 

1.2 

Noise stationarity 

16.9 

17.0 

1.4 

FPU Temperature 

16.7 

16.7 

0.84 

Telescope pointing 

15.7 

15.8 

3.0 

Passing data 

15.4 

15.6 

0.55 


Bicep 2 Results and Instrument Papers the map depth for the 
Keck Array 2012-2013 150 GHz data is 74 nK in nominal 
square-degree pixels (4.4 /iKarcmin) over an effective area of 
390 square degrees for a total sensitivity of 2.6 nK. An equiv¬ 
alent way of expressing the sensitivity of the data set is the 
survey weight W = \/s^ =150,000 where s is the total 
sensitivity. This expression is useful because it scales linearly 
with integration time, number of detectors, and statistical sen¬ 
sitivity to r. 

5. LOW LEVEL DATA REDUCTION, MAP MAKING AND 
SIMULATIONS 

5.1. Analysis pipeline 

The Keck Array and Bicep2 data analysis uses the same 
code and proceeds in parallel, enabling cross checks between 
these independent data sets. The process used here is exactly 
the same as described in the Bicep2 Results Paper—a sum¬ 
mary follows. 

5.2. Low level reduction 

As detailed in the Bicep2 Results Paper, the first step in the 
low level reduction is to deconvolve the temporal response of 
the instrument. The TES detectors have a fast response of 
~1 ms that can be ignored. Both the MCE and the gcp apply 
low pass filters to the data which must be accounted for in the 
deconvolution. 

At this stage, relative calibration is accomplished by divid¬ 
ing the time streams by the individual detector gains derived 
from elevation nods. The data are then multiplied by the me¬ 
dian gain across the array in order to remove dependence on 
atmospheric variation. 

5.3. Pairmaps 

The sum and difference of each detector pair is taken. Each 
half-scan is subjected to third order polynomial filtering to 
remove atmospheric variation. In addition, the mean of the 
half-scans over the scanset is subtracted to remove any scan- 
synchronous contamination. The pointing of each detector 
pair is reconstructed from the telescope pointing model and 
per-pair offset angles refined by regressing per-channel maps 
against a WMAP5 template. At this point the time stream 
data of each pair is binned into a rectangular grid of pixels 
forming per scanset “pairmaps”. We also sample and bin the 
Planck 143 GHz temperature map and its derivatives to use in 


the deprojection of beam systematics. Eor further details see 
the Bicep2 Results Paper and the Systematics Paper. 

5.4. Full maps 

Einally, the pairmaps are coadded into final full maps, and 
also into various pairs of jackknife splits. The deprojection 
templates are fit and removed during this process. Absolute 
calibration of the maps is performed by comparing the power 
spectrum of the temperature map with the Planck 143 GHz 
map as described in the Instrument Paper. 

Eigure 3 shows the resulting temperature and Stokes Q and 
U polarization maps for the Keck Array 150 GHz data from 
the 2012-2013 observing seasons. The left column shows 
the final maps and the right hand column shows a difference 
(jackknife) map which is consistent with noise alone. The ver¬ 
tical/horizontal stripes in the Q maps and the diagonal stripes 
in the U maps are characteristic of the E-mode polarization 
signal, which dominates the maps. 

5.5. Simulations 

We create signal and noise simulations exactly as described 
in the Bicep 2 Results Paper. We generate realizations of 
noise by randomly flipping the signs of the pairmaps when 
co-adding to full maps. Several kinds of signal simulations 
are made by resampling input maps from the synf ast pro¬ 
gram (part of the Healpix^^ package (Gorski et al. 2005)). 
The simulated data are then binned into pairmaps and com¬ 
bined to full maps exactly in parallel with the treatment of the 
real data. 

6. RESULTS 
6.1. Power spectra 

The maps are converted into angular power spectra exactly 
as described in the Bicep2 Results Paper. The matrix based 
purification of the Q and U maps is performed prior to inver¬ 
sion to form Z?-modes to avoid E to B mixing due to the sky- 
cut and filtering. A variant of the MASTER procedure (Hivon 
et al. 2002) is used to noise debias the auto spectra and correct 
for the effects of the time stream filtering. 

The resulting power spectra for the Keck Array using the 
150 GHz data from the 2012-2013 observing seasons is 
shown in Eigure 4, along with a temporal jackknife. The B- 
mode power spectrum is inconsistent with ACDM cosmology 

http: //healpix. sourcef orge . net/ 
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Fig. 2.— The Keck Array 2012-2013 150 GHz data set. The top panel shows the fraction of calendar time the telescopes were observing. The red line 
represents the fraction of time preserved after all data selection cuts are implemented. The middle panel shows the instantaneous sensitivity of the full Keck 
Array. The bottom panel shows the cumulative map depth as calculated over each phase (10 hr of data). 


without foregrounds at 5 .Oct (over the first 5 band powers). 
Although the noise is lower than for Bicep 2, the first two 
band power values are also lower so the significance is some¬ 
what smaller. 

An overall rotation is applied to the maps to minimize the 
high-£ TB and EB spectra (Kaufman et al. 2014). For the 
Keck Array 2012H-2013 data this adjustment is « -0.5°. This 
rotation makes no practical difference to the B-mode power 
spectrum. 

6.2. E and B maps 

F-mode and B-mode maps can be made by performing an 
inverse Fourier transform as shown in Figure 5. The maps cre¬ 
ated are inherently apodized, as the B-mode and B-mode com¬ 
ponents are generated from apodized Q and U maps. These 
are compared to a lensed-ACDM-hnoise simulation. 

6.3. Internal consistency tests 

The Keck Array data was split in 16 different ways to test for 
internal consistency. The motivations behind these splits are 
described in the Bicep2 Results Paper and the Systematics 
Paper. If a contaminating signal exists in only one half of the 
data split, then it should show up with as much significance in 
the jackknife as in the signal map. However, some jackknives 
are more sensitive to certain systematics than the signal map 


because of inherent cancellation effects which operate in the 
full map. Each of the jackknife categories is summarized be¬ 
low. 

The first set of jackknives probes for systematics which dif¬ 
fer between different subsets of channels. This includes divi¬ 
sion in the multiplexing system, as well as divisions in the fo¬ 
cal plane layout: tile, focal plane inner/outer, tile inner/outer, 
mux row, mux column. As is documented in the Beams Paper, 
there are systematics that are highly dependent on the position 
of the detector in the focal plane. For instance, the ellipticity 
of the beams is greater in the detectors near the outside of the 
focal plane than the inside. 

The next set of jackknives is temporal. This includes both 
the longest time scale of 2012 data versus 2013 data, and the 
shortest timescale of left going scans versus right going scans. 
Owing to the changes between the 2012 and 2013 observing 
seasons, an early/late season jackknife acts as an alternative 
temporal split. The first is sensitive to the effects of different 
observing schedules and detectors changed between seasons, 
while the second is sensitive to detector time constants. 

Another set of jackknives is based on external contamina¬ 
tion. This includes the azimuth jack, which divides the data 
based on the direction the telescope is pointed with respect 
to the ground (see Table 2). One half of this jackknife is data 
taken in the direction of the main South Pole Station and asso- 
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Fig. 3.— Keck Array T, Q, U maps. The left column shows the basic signal maps with 0.25° pixelization as output by the reduction pipeline. The right column 
shows difference (jaclcknife) maps made with the first halves of the 2012 and 2013 seasons and the second halves. No additional filtering other than that imposed 
by the instrument beam (FWHM 0.5°) has been done. Note that the structure seen in the Q and U signal maps is as expected for an iJ-mode dominated sky. 


elated operations, while the other half points into the desolate 
Antarctic plateau. 

A set of jackknives that particularly amplifies the differen¬ 
tial beam properties is the deck rotation jackknives. As is 
described in the Bicep2 Systematics Paper, a 180° deck ro¬ 
tation cancels out differential pointing. The deck jackknife, 
which differences the 180° rotations, amplifies the leakage by 
an order of magnitude in comparison to that present in the 
fully coadded data. The Keck Array also started taking data at 
90° compliment deck rotations in the 2013 observing season, 
and this jackknife is sensitive to differential gain or differen¬ 
tial beam width leakage. The alternative deck jackknife is 
defined to be the difference of the 90° rotations for 2013. In 
this special case, the statistics for the 2012 and 2013 data are 
separate. 

Maps are made from each half of the data split and then dif¬ 
ferenced. The differenced maps are divided by a factor of two 
in order to keep the noise amplitudes equivalent to the signal 
map. The consistency with lensed-ACDM-i-noise simulations 
is calculated with a simple statistic; 

X^= (d-(m))'^D“'(d-(m)) (6) 

where d is the vector of observed band power values, (m) 
is the mean of the lensed-ACDM-i-noise simulations (except 
where alternative signal models are considered), and D is the 
band power covariance matrix as evaluated from those simu¬ 


lations. 

A X statistic is also considered to probe for sets of band 
powers which are systematically above or below the expecta¬ 
tion. This is defined as; 

where the di are the observed band power values and (m,) 
and (Tm; are the mean and standard deviation of the lensed- 
ACDM-i-noise simulations. 

For each of these statistics, we calculate the probability to 
exceed (PTE) the observed value by comparing to the values 
obtained in the 500 lensed-ACDM-i-noise simulations. Since 
the distribution of the band powers of the auto spectra is ap¬ 
proximately x^ distribution, there is some non-Gaussianity to 
the statistics. In particular, the lowest band power only has 
9 effective modes which will increase the tails of the distri¬ 
bution. However, by calculating the PTE against the simula¬ 
tions, any non-Gaussianity is fully reflected in the PTE value. 
The PTE for the x and using band powers 1-5 and 1-9 is 
given in Table 4. Note that these statistics are correlated (es¬ 
pecially along each row of the table). The distribution of the 
PTE values is shown in Eigure 6. 

7. SYSTEMATICS 
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Fig. 4.— Keck Array power spectrum results for signal (black points) and early/late season jackknife (blue points). The solid red curves show the lensed- 
ACDM theory expectations.The eiTor bars ai'e the standard deviations of the lensed-ACDM+noise simulations and hence contain no sample variance on any 
additional signal component. The probability to exceed (PTE) the observed value of a simple statistic for the 9 band powers is given (as evaluated against 
the simulations). Note that the band powers of the auto spectra of the simulations are approximately distributed, with the lowest i-bin only containing ~9 
effective degrees of freedom. This increases the probability of an outlier point in comparison to a Gaussian distribution. The observed distribution for the cross 
spectra are more symmetric than the x^ distributions but have similarly increased tails. This is fully reflected in the quoted PTE value. Also note the very different 
y-axis scales for the jackknife spectra (other than BB). See the text for additional discussion of the BB spectrum. (Note that the calibration procedure uses EB to 
set the overall polarization angle so TB and EB as plotted above cannot be used to measure astrophysical polarization rotation.) 


Experimental systematics can create false B-mode polariza¬ 
tion and must be shown to be tightly controlled. The system¬ 
atics in Bicep 2 were fully explored in the Systematics Paper 
and were shown to be below the level equivalent to r=0.003- 
0.006. These limits were derived from forward simulations 
of the measured instrumental properties. If a given property 
did not have a measured level, appropriate upper limits were 
used. 

The beam systematics in the Keck Array are expected to be 
below those of Bicep2 because of the larger number of de¬ 
tectors (increased averaging down of incoherent effects), and 
the increased number of receiver orientations (both instanta¬ 
neously due to the “clocking” of the five receivers at 72° in¬ 
crements around the boresight, and, in 2013, the increased 
number of deck angles of observation). 

As described in the Bicep2 Results and Systematics Papers 
we produce simulated time streams by convolving an input 
temperature map with high precision per channel measure¬ 
ments of the actual beam shapes. We then pass these simu¬ 
lated time streams through the mapping process, including all 
filtering and deprojection, to assess the residual contamina¬ 
tion due to beam non-ideality. The results are shown in Fig¬ 
ure 7. The beam maps for the Keck Array do not provide as 
uniform and redundant coverage of all detectors as those for 
Bicep2 and additional analysis is required to construct com¬ 


posite beam maps that have consistently high signal-to-noise 
and are free of artifacts from the beam mapping measurement. 
For the purposes of the current paper we use the preliminary 
beam map results to set an upper limit on the residual contam¬ 
ination, as indicated by the down arrows in the figure. 

Other forms of systematic contamination were considered, 
such as electromagnetic interference (EMI) contamination, 
magnetic pickup, thermal pickup, and detector pointing. Each 
of these systematics was quantified to be below the Bicep 2 
level when averaged over the entire array, and thus safely ig- 
norable. 

8. CONSISTENCY WITH Bicep2 

The resulting B-mode power spectra of Bicep2 and the 
Keck Array can be compared to assess the compatibility of 
the two sets of results. Although there is much that the two 
experiments share in terms of hardware, design, and location, 
there are also potential systematics that are different; the bulk 
refrigeration system, the ground shield, and the time at which 
the observations occurred. Comparing the results is a power¬ 
ful additional systematics check. The power spectra for both 
the Keck Array, Bicep2 and the cross between the two are 
shown in the upper panel of Figure 8. 

A rigorous comparison can be done in two ways: directly 
comparing the maps and comparing the auto and cross power 
spectra. The latter can be a more powerful comparison if the 
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Keck Array : E signal Simulation: E from lensed-ACDM+noise 



Right ascension [deg.] 


Fig. 5.— Left: Keck Array apodized iJ-mode and B-mode maps filtered to 50 < ^ < 120. Right: The equivalent maps for the first of the lensed-ACDM+noise 
simulations. The color scale displays the F-mode scalar and B-mode pseudoscalar patterns while the lines display the equivalent magnitude and orientation of 
linear polarization. Note that the B-mode and B-mode maps use different color/length scales. 
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Fig. 6.— Distributions of the jackknife ^nd x PTE values over the tests 
and spectra given in Table 4. 


maps have different noise levels—since Bicep 2 and the Keck 
Array have comparable noise levels, all four methods (map 
and the three combinations of auto and cross spectra) have 
approximately equal statistical power. 

To test the compatibility of the resulting band powers with 
null we compare them to the differences of signal+noise sim¬ 
ulations which share common input skies. In such tests it is 
necessary that the simulations contain power roughly equal 



Multipole 


Fig. 7.— BB spectra from T-only input simulations using the measured 
per channel beam shapes compared to the lensed-ACDM+r = 0.2 spectnim. 
From top to bottom the curves are (i) no deprojection, (ii) deprojection of dif¬ 
ferential pointing only (<5jc+(5y), (iii) deprojection of differential pointing and 
differential gain of the detector pairs ((5.x' + 5v + 5g), (iv) adding deprojection 
of differential beam width ((5x + (5y + (5g + (5(j), and (v) differential pointing, 
differential gain, and differential ellipticity ((5x+5y + (5g + (5/? + (5c‘). This last 
curve represents an upper limit only to the residual contamination. 

to the real sky as the cross terms between signal and noise 
increase the fluctuation of the differences even for perfectly 
common sky coverage. For example, the (un-debiased) auto 
spectrum of a map M composed of a signal S and noise N can 
be written 

MxM=(S-\-N)x(S-\-N) = SxS-\-2(SxN)-\-NxN. (8) 
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TABLE 4 

Jackknife PTE values from and x (sum of deviation) tests 


Jackknife 

Band powers 

1-5 x" 

Band powers 

1-9 x' 

Band powers 

1-5 X 

Band powers 

1-9 X 

Deck Jackknife 




EE 

0.613 

0.924 

0.898 

0.776 

BB 

0.743 

0.880 

0.685 

0.375 

EB 

0.820 

0.986 

0.309 

0.429 

Scan Dir jackknife 




EE 

0.561 

0.415 

0.788 

0.898 

BB 

0.924 

0.691 

0.601 

0.180 

EB 

0.168 

0.453 

0.938 

0.886 

Early/Late Season jackknife 




EE 

0.287 

0.255 

0.896 

0.998 

BB 

0.982 

0.960 

0.621 

0.796 

EB 

0.711 

0.667 

0.170 

0.609 

Year Split jackknife 




EE 

0.343 

0.641 

0.918 

0.956 

BB 

0.856 

0.940 

0.695 

0.657 

EB 

0.747 

0.547 

0.353 

0.798 

Tile jackknife 




EE 

0.042 

O.IIO 

0.431 

0.782 

BB 

0.573 

0.715 

0.118 

0.499 

EB 

0.451 

0.691 

0.940 

0.932 

Phase jackknife 




EE 

0.826 

0.824 

0.743 

0.309 

BB 

0.036 

0.184 

0.489 

0.343 

EB 

0.026 

0.058 

0.980 

0.914 

Mux Col jackknife 




EE 

0.804 

0.142 

0.543 

0.080 

BB 

0.471 

0.760 

0.291 

0.295 

EB 

0.144 

0.206 

0.585 

0.840 

Alt Deck jackknife 2012 




EE 

0.673 

0.884 

0.641 

0.377 

BB 

0.579 

0.685 

0.569 

0.784 

EB 

0.152 

0.112 

0.389 

0.609 

Alt Deck jackknife 2013 




EE 

0.489 

0.583 

0.269 

0.507 

BB 

0.583 

0.822 

0.794 

0.826 

EB 

0.549 

0.441 

0.577 

0.575 

Mux Row jackknife 




EE 

0.942 

0.645 

0.421 

0.776 

BB 

0.363 

0.371 

0.625 

0.786 

EB 

0.214 

0.345 

0.030 

0.112 

Tile/Deck jackknife 




EE 

0.204 

0.052 

0.910 

0.240 

BB 

0.487 

0.745 

0.986 

0.830 

EB 

0.491 

0.351 

0.146 

0.603 

Focal Plane inner/outer jackknife 



EE 

0.253 

0.475 

0.108 

0.064 

BB 

0.637 

0.587 

0.074 

0.158 

EB 

0.044 

0.226 

0.996 

0.994 

Tile top/bottom jackknife 




EE 

0.573 

0.397 

0.924 

0.910 

BB 

0.289 

0.194 

0.198 

0.782 

EB 

0.172 

0.353 

0.884 

0.954 

Tile inner/outer jackknife 




EE 

0.707 

0.663 

0.399 

0.387 

BB 

0.303 

0.663 

0.719 

0.786 

EB 

0.958 

0.655 

0.315 

0.102 

Moon jackknife 




EE 

0.192 

0.433 

0.395 

0.385 

BB 

1.000 

0.387 

0.339 

0.505 

EB 

0.667 

0.705 

0.794 

0.289 

A/B offset best/worst 




EE 

0.483 

0.521 

0.804 

0.926 

BB 

0.443 

0.367 

0.042 

0.407 

EB 

0.497 

0.677 

0.489 

0.397 



Fig. 8.— Upper: The Keck Array BB auto spectrum, the Bicep2 auto spec¬ 
trum, and the cross spectrum taken between the two. The inner error bars 
are the standard deviation of the lensed-ACDM+noise simulations, while the 
outer en'or bars also contain excess power at low-£. (For clarity the Keck 
Array and cross spectrum points are offset horizontally.) Lower: Four com¬ 
patibility tests between the B-modes measured by Bicep2 and Keck Array. 
The “map jack” takes the difference of the Q and U maps, divides by a factor 
of two, and calculates the BB spectrum. The other three sets of points are 
the differences of the spectra shown in the upper panel divided by a factor of 
four. In each case the eiTor bai's are the standai'd deviation of the pairwise dif¬ 
ferences of signal+noise simulations which share common input skies. Com¬ 
parison of any one of these sets of points with null is an appropriate test of the 
compatibility of the experiments—see text for details. All tests show good 
consistency between Bicep2 and Keck Array, particularly in the lowest five 
bandpowers. 

The difference of such auto spectra between experiments with 
common sky coverage is then 

MixMi -M 2 XM 2 = 2(SxNi-SxN2)+NixNi-N 2 XN 2 (9) 

where Mi and M 2 refer to the first and second experiment. 
The signal auto spectrum SxS cancels out. However, the cross 
terms between the signal and noise 2(SxNi -SXN 2 ) do not 
cancel, and they increase the fluctuations between the two ex¬ 
periments over the noise-only case in proportion to the com¬ 
mon signal. To account for this extra variance, we use signal 
simulations with additional power that matches the amplitude 
of the observed signal in excess of ACDM in band powers 1- 
5. (The origin of the extra signal over ACDM is not important 
here—only its approximate amplitude.) The results are shown 
in Figure 8. 

We then proceed to calculate the PTE of the x and statis¬ 
tics versus the simulated distributions using the same spectra 
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TABLE 5 

BicnvnKeck Array COMPATIBILITY TEST PTE VALUES FROM and X 
(SUM OF DEVIATION) TESTS 


Jackknife 

Band powers 

1-5 x" 

Band powers 

1-9 x" 

Band powers 

1-5 X 

Band powers 

1-9 X 

Map jackknife 




EE 

0.034 

0.048 

0.106 

0.028 

BB 

0.561 

0.695 

0.054 

0.018 

EB 

0.741 

0.754 

0.405 

0.651 

Spectral jackknife B2-cross 




EE 

0.112 

0.092 

0.068 

0.078 

BB 

0.687 

0.387 

0.052 

0.008 

EB 

0.555 

0.224 

0.212 

0.234 

Spectral jackknife B2-Keck 




EE 

0.138 

0.128 

0.066 

0.126 

BB 

0.920 

0.485 

0.200 

0.022 

EB 

0.511 

0.214 

0.210 

0.200 

Spectral jackknife cross-Keck 




EE 

0.176 

0.204 

0.074 

0.202 

BB 

0.880 

0.966 

0.643 

0.435 

EB 

0.361 

0.437 

0.443 

0.188 

and band 

power ranges as in §6.3, and give the results in Ta- 


ble 5. In both the figure and the table we note the effect of 
the two band powers at f « 220 that are high with respect to 
lensed-ACDM in B2xB2 (as noted in the Bicep2 Results Pa¬ 
per) but not in KxK and B2xK—as expected these also show 
up in the map difference. Again note that the PTE values are 
correlated (both along and between rows of the table) so over¬ 
interpretation should be avoided. Our conclusion is that the 
Bicep2 and Keck Array data are consistent—especially in the 
lowest five band powers where an IGW contribution would be 
strongest. 

9. COMBINATION WITH Bicep2 

Having shown that the Keck Array results are consistent 
with Bicep2 we now proceed to combine the maps by adding 
the accumulation quantities (equivalent to a noise weighted 
combination of the maps). This results in Q and U maps 
which have a depth of 57 nKdeg (3.4 pKarcmin) over an 
effective area of 400 square degrees. Following §4.3, the map 
depth and effective area are combined for a total sensitivity of 
2.0 nK and a total survey weight of 250,000 

The observation regions and strategies are sufficiently sim¬ 
ilar that it is found empirically using simulations that the pu¬ 
rification matrix of either experiment delivers adequate B- 
mode purity when applied to the combined map (with con¬ 
tamination equivalent to r < 10“^). 

The final BB spectrum is shown in Figure 9 and is incon¬ 
sistent with the lensed-ACDM expectation at > 6cr (for ei¬ 
ther band powers 1-5 or 1-9). The lensed-ACDMn-noise er¬ 
ror bars as plotted are approximately a factor two smaller 
than those of the previous Bicep2 only results—saturation 
on the (small) sample variance of the lensing component is 
occurring—the noise component is a factor 2.3 times smaller. 
All the spectra (including TT, EE etc.) are available for 
download at http : / /bicepkeck . org/ together with the 
ancillary data, noise information etc., required to use them. 

10. CONCLUSIONS 

We have presented the Keck Array instrument and the 
2012-2013 (150 GHz) data set. The instantaneous instrumen¬ 
tal sensitivity of 9.5 is the best reported to date. 

The same area of sky as previously observed by Bicep2 was 



Fig. 9.— The BB power spectrum of combined Bicep2 and Keck Ar¬ 
ray maps. The inner error bai's are the standard deviation of the lensed- 
ACDM+noise simulations, while the outer error bars also contain excess 
power at low-^. 

mapped to a depth in Q and U of 74 nKdeg (4.4 /iKarcmin). 
The resulting Keck Array power spectra are consistent with 
lensed-ACDM except for an excess at degree angular scales in 
BB which has a significance of 5.0cr. Extensive jackknife tests 
argue against a systematic origin for the signal, and further 
statistical tests indicate that the maps and spectra are consis¬ 
tent with the previous Bicep2 results. Finally the two sets of 
maps are combined to produce maps with noise of 57 nKdeg 
(3.4 pKarcmin) over an effective area of 400 deg^ for a sur¬ 
vey weight of 250,000 The final BB spectrum is in¬ 

consistent with lensed-ACDM at a significance of > 6a. The 
combined map results (for all spectra) are available for down¬ 
load. There does not appear to be any reason to consider the 
Bicep2 results as more reliable than the Keck Array results or 
vice versa. We therefore emphasize that we regard the com¬ 
bined results as the best available data set at this time. 

The origin of the excess power shown in Figure 9, and pre¬ 
viously reported in the Bicep 2 Results Paper, has been ex¬ 
tensively debated in the literature (Mortonson & Seljak 2014; 
Flauger et al. 2014; Fuskeland et al. 2014). Recently, con¬ 
crete information on the strength of polarized dust emission 
at high galactic latitude has become available in Planck Col¬ 
laboration Int. XXX (2014). It appears that dust emission is 
a significant contribution to the signal observed by Bicep 2 
and the Keck Array. Therefore, in an upcoming paper, the 
Bicep 2 and combined maps are cross correlated with Planck 
maps of the same region to constrain the dust contribution to 
the observed signal. 

During the 2014 season two of the Keck Array receivers 
operated at 95 GHz and a future analysis will use this data 
to further constrain the dust contribution. In the 2015 sea¬ 
son Bicep3 will provide increased sensitivity at 95 GHz and 
operation of Keck Array receivers at 220 GHz is also planned. 
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in particular the South Pole Station without whose help this 
research would not have been possible. Most special thanks 
go to our heroic winter-overs Robert Schwarz and Steffen 
Richter. We thank all those who have contributed past efforts 
to the ^lCE¥~Keck Array series of experiments, including the 
Bicep 1 team. 


REFERENCES 


Ade, P. A. R., Pisano, G., Tucker, C., & Weaver, S. 2006, Proc. SPIE, 6275, 
62750U 

Aikin, R. W., Ade, P. A., Benton, S., et al. 2010, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7741, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 0 
Barkats, D., Bischoff, C., Farese, P, et al. 2005, Astrophys. J. Lett., 619, 
L127 

Battistelli, E. S., et al. 2008, J. Low Temp. Phys., 151, 908 
Bennett, C. L., Larson, D., Weiland, J. L., et al. 2013, Astrophys. J. Suppl., 
208, 20 

Benton, S. J., Ade, P. A., Amiri, M., et al. 2014, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 9145, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 0 
BicepI Collaboration, Barkats, D., Aikin, R., et al. 2014, The Astrophysical 
Journal, 783, 67 

Bicep2 Collaboration I. 2014, Phys. Rev. Lett., 112, 241101 
Bicep2 Collaboration 11. 2014, The Astrophysical Journal, 792, 62 
Bicep2 Collaboration III. 2015, In prep. 

Bicep2, Keck Array, and SPIDER Collaborations. 2015, In prep. 

Bicep2 / Keck Array Collaborations IV. 2015, In prep. 

Bischoff, C., Hyatt, L., McMahon, J. J., et al. 2008, Astrophys. J., 684, 771 
Brown, M. L., Ade, R, Bock, J., et al. 2009, Astrophys. J., 705, 978 
Caligiuri, J., & Kosowsky, A. 2014, Physical Review Letters, 112, 191302 
Chiang, H. C., Ade, P. A. R., Barkats, D., et al. 2010, Astrophys. J., 711, 
1123 

Crites, A. T., Henning, J. W., Ade, P. A. R., et al. 2014, ArXiv e-prints, 
arXiv:1411.1042 

de Korte, P. A. J., Beyer, J., Deiker, S., et al. 2003, Review of Scientific 
Instruments, 74, 3807 

Dodelson, S. 2014, Physical Review Letters, 112, 191301 
Duband, L., & Collaudin, B. 1999, Cryogenics, 39, 659 
Flauger, R., Hill, J. C., & Spergel, D. N. 2014, arXiv:1405.7351vl 
Fuskeland, U., Wehus, 1. K., Eriksen, H. K., & Nsss, S. K. 2014, Astrophys. 
J., 790, 104 

Gain, S., Benabed, K., Bouchet, R, et al. 2014, Phys. Rev. D, 90, 063504 
Gorski, K. M., Hivon, E., Banday, A. J., et al. 2005, Astrophys. J., 622, 759 
Hanson, D., Hoover, S., Crites, A., et al. 2013, Physical Review Letters, 111, 
141301 

Hinderks, J. R., Ade, R, Bock, J., et al. 2009, Astrophys. J., 692, 1221 
Hivon, E., Gorski, K. M., Netterfield, C. B., et al. 2002, Astrophys. J., 567, 2 
Irwin, K., & Hilton, G. 2005, in Topics in Applied Rhysics, Vol. 99, 
Cryogenic Rarticle Detection, ed. C. Enss (Springer Berlin Heidelberg), 
63-150 

Kamionkowski, M., Kosowsky, A., & Stebbins, A. 1997, Rhysical Review 
Letters, 78, 2058 


Karkare, K. S., Ade, R. A. R., Ahmed, Z., et al. 2014, in Society of 

Rhoto-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 
9153, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, 3 

Kaufman, J. R, Miller, N. J., Shimon, M., et al. 2014, Rhys. Rev. D, 89, 
062006 

Kernasovskiy, S., et al. 2012, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 8452, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series 
Kovac, J. M., Leitch, E. M., Pryke, C., et al. 2002, Nature, 420, 772 
Leitch, E. M., Kovac, J. M., Pryke, C., et al. 2002, Nature, 420, 763 
Montroy, T. E., Ade, P. A. R., Bock, J. J., et al. 2006, Astrophys. J., 647, 813 
Mortonson, M. J., & Seljak, U. 2014, arXiv: 1405.5857 
Naess, S., Hasselfield, M., McMahon, J., et al. 2014, JCAP , 10, 7 
Page, L., Hinshaw, G., Komatsu, E., et al. 2007, Astrophys. J. Suppl., 170, 
335 

Planck Collaboration Int. XXX. 2014, Astr. & Astroph., arXiv: 1409.5738 
Planck Collaboration results 1. 2014, Astr. & Astroph., 571, Al 
POLARBEAR Collaboration, Ade, P. A. R., Akiba, Y., et al. 2014a, 
Astrophys. J., 794, 171 

—. 2014b, Physical Review Letters, 112, 131302 
—. 2014c, Physical Review Letters, 113, 021301 
Polnarev, A. G. 1985, Sov. Ast., 29, 607 
Pryke, C., Ade, R, Bock, J., et al. 2009, Astrophys. J., 692, 1247 
Quiet Collaboration, Bischoff, C., Brizius, A., et al. 2011, Astrophys. J., 
741, 111 

Quiet Collaboration, Araujo, D., Bischoff, C., et al. 2012, Astrophys. J., 
760, 145 

Readhead, A. C. S., Myers, S. T., Pearson, T. J., et al. 2004, Science, 306, 
836 

Rocha, G., Trotta, R., Martins, C. J. A. R, et al. 2004, Mon. Not. Roy. 

Astron. Soc., 352, 20 
Seljak, U. 1997, Astrophys. J., 482, 6 

Seljak, U., & Zaldarriaga, M. 1997, Physical Review Letters, 78, 2054 
Sheehy, C. D., et al. 2010, Millimeter, Submillimeter, and Far-Infrared 
Detectors and Instrumentation for Astronomy V, 7741, 7741IG 
Sievers, J. L., Achermann, C., Bond, J. R., et al. 2007, Astrophys. J., 660, 
976 

Stiehl, G., Cho, H.-M., Hilton, G., et al. 2011, Applied Superconductivity, 
IEEE Transactions on, 21, 298 

Story, K., Leitch, E., Ade, R, et al. 2012, Proc. SPIE Int. Soc. Opt. Eng., 
8451, 845 lOT 

Takahashi, Y. D., Ade, P. A. R., Barkats, D., et al. 2010, Astrophys. J., 711, 
1141 

Wu, J. H. R, Zuntz, J., Abroe, M. E., et al. 2007, Astrophys. J., 665, 55 
Zaldarriaga, M., & Seljak, U. 1998, Phys. Rev. D, 58, 023003 



